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Abstract 

From the 2002 data taking with a neutral kaon beam extracted from the CERN-SPS, the 
NA48/1 experiment observed 97 S° — > E + ^~77 p candidates with a background contami- 
nation of 30.8 ± 4.2 events. From this sample, the BR(S° — » E + yU~77 M ) is measured to be 
(2.17±0.32 stat ±0.17 syst ) x 10- 6 . 
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1 Introduction 

The study of hadron beta decays gives important information on the interplay between the 
weak interaction and hadron structure determined by the strong interaction. In this framework, 
measurements on H° semileptonic decays and on the related parameters are fundamental to 
further increase our knowledge on the constituents of the baryon octet. In particular, a clear 
evidence for the decay H° — > and a measurement of its branching ratio will add one 

more constraint to the theoretical frameworks [2J O HI EJ E] built to explain the behaviour of the 
baryon semileptonic decays. 

In the present article the branching ratio of the semileptonic decay H° — > S + [i^v^ is 
measured by normalizing to the analogue decay with an electron in the final state, already 
studied by the NA48/1 collaboration [7J. The similar topologies of the final states of the two 
semileptonic decays allowed the same trigger conditions to be used for both data samples. The 
selection criteria are also similar between the two channels and only differ for the identification 
of the charged lepton and for cuts related to background rejection. This decay had already 
been observed by the KTeV collaboration [8], with a sample of 9 events and a branching ratio 
measurement of (4.7^'g) x 10~ 6 . 

2 Beam 

The experiment was performed in 2002 at the CERN SPS accelerator and used a 400 GeV 
proton beam impinging on a Be target to produce a neutral beam. The spill length was 4.8 s 
out of a 16.2 s cycle time. The proton intensity was fairly constant during the spill with a mean 
of 5 x 10 10 particles per pulse. 

For this measurement, only the K$ target station of the NA48 double Ks/Kl beam line [S] 
was used to produce the neutral beam. In this configuration, the Kl beam was blocked and an 
additional sweeping magnet was installed to deflect charged particles away from the defining 
section of the Kg collimators. To reduce the number of photons in the neutral beam originating 
primarily from 7r° decays, a 24 mm thick platinum absorber was placed in the beam between 
the target and the collimator. A pair of coaxial collimators, having a total thickness of 5.1 m, 
the axis of which formed an angle of 4.2 mrad to the proton beam direction, selected a beam of 
neutral long-lived particles (K$, Kl, A , H°, n and 7). The target position and the production 
angle were chosen in such a way that the beam axis was hitting the center of the electromagnetic 
calorimeter. 

In order to minimize the interaction of the neutral beam with air, the collimator was 
immediately followed by a 90 m long evacuated tank terminated by a 0.3% Xq thick Kevlar 
window. The NA48 detector was located downstream of this region. 

On average, about 1.4 x 10 4 H° per spill, with an energy between 70 and 220 GeV, decayed 
in the fiducial decay volume. 

3 Detector 

The detector was designed for the measurement of Re(e' /e), and a detailed description of 
the experimental layout is available at [9]. In the following sections a short description of the 
main detectors is reported. 

3.1 Tracking 

The detector included a spectrometer housed in a helium gas volume with two drift cham- 
bers before and two after a dipole magnet with an horizontal transverse momentum kick of 265 
MeV/c. Each chamber had four views (x, y, u, v), each of which had two sense wire planes. 
The resulting space points were typically reconstructed with a resolution of ~ 150 /im in each 
projection. The spectrometer momentum resolution is parameterized as: 

a p /p = 0.48% 0.015% x p 
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where p is in GeV/ c. This gave a resolution of 3 MeV/ c 2 when reconstructing the kaon mass in 
K° — > tt + it~ decays. The track time resolution was ~ 1.4 ns. 

3.2 Electromagnetic Calorimetry 

The detection and measurement of the electromagnetic showers were achieved with a liquid 
krypton calorimeter (LKr), 27 radiation lengths deep, with a ~2 cm x 2 cm cell cross-section. 
The energy resolution, expressing E in GeV, is parameterized as [9]: 

a(E)/E = 3.2%/\/£0 9%/£ 0.42% 

The transverse position resolution for a single photon of energy larger than 20 GeV was 
better than 1.3 mm, and the corresponding mass resolution for the reconstructed ir° mass (77 
decay) was ~1 MeV/c 2 . The time resolution of the calorimeter for a single shower was better 
than ~ 300 ps. 

3.3 Scintillator Detectors and Muon Detector 

A scintillator hodoscope (CHOD) was located between the spectrometer and the calorime- 
ter. It consisted of two planes, segmented in horizontal and vertical strips and arranged in 
four quadrants. Further downstream there was an iron-scintillator sandwich hadron calorimeter 
(HAC), followed by muon counters consisting of three planes of scintillator, each shielded by an 
80 cm thick iron wall. The first two planes MIX and M1Y were the main muon counters and 
had 25 cm wide horizontal and vertical scintillator strips respectively, with a length of 2.7 m. 
The third plane, M2X, had horizontal strips 44.6 cm wide, and was mainly used to measure the 
efficiency of the MIX and M1Y counters. The central strip in each plane was divided into two 
sections separated by a gap of 21 cm, in order to accommodate the beam pipe. The fiducial vol- 
ume of the experiment was principally determined by the LKr calorimeter acceptance, together 
with seven rings of scintillation counters (AKL) used to veto activity outside this region. 

4 Trigger 

The trigger system used for the on-line selection of H° semileptonic decays mainly consisted 
of two levels of logic. Level 1 (LI) was based on logic combinations of fast signals coming from 
various sub-detectors. It required hits in the CHOD and in the first drift chamber compatible 
with at least one and two tracks respectively, no hit in the AKL veto system and a minimum 
energy deposition in the calorimeters. This last requirement was 15 GeV for the energy re- 
constructed in the LKr calorimeter or 30 GeV for the summed energy in the electromagnetic 
and hadronic calorimeters. The output rate of the LI stage was about 50 kHz. The average LI 
efficiency, measured with H° — > Air events of energy greater than 70 GeV, was found to be 
98.65 ± 0.03%. 

Level 2 (L2) consisted of a set of 300 MHz processors that reconstructed tracks and ver- 
tices from hits in the drift chambers and computed relevant physical quantities. The L2 trigger 
required at least two tracks with a closest distance of approach of less than 8 cm in space and 
a transverse separation greater than 5 cm in the first drift chamber. Since the signature of the 
H° /3-decay involves the detection of an energetic proton from the subsequent S + — > pn decay, 
the ratio between the higher and the lower of the two track momenta was required to be larger 
than 3.5. Rejection of the overwhelming A — > pn~ and K$ —> ir + ir~ decays was achieved by 
applying stringent invariant mass cuts against these decays. The output L2 trigger rate was 
about 2.5 kHz. The efficiency of the L2 trigger stage with respect to Level 1, averaged over the 
2002 run, was measured to be (83.7±2.2)% for H° /3-decays, mainly limited by wire inefficiencies 
in the drift chambers. 
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5 Offline selection 

The identification of the H° — > H + fj,~Uu channel was performed using the subsequent decay 
S + —?• p7T° with 7T° — > 77. The final state consists of a proton and a muon, giving two tracks in 
the spectrometer, two photons producing clusters in the LKr calorimeter and one unobserved 
anti-neutrino. The decay S° — > Tj + 1~V^ is the only source of S + particles in the neutral beam 
since the two-body decay H° — > £ + 7r~ is kinematically forbidden. Thus, the signal events were 
identified by requiring an invariant pir° mass consistent with the nominal S + mass value. 

The S + decay was reconstructed using a positive charged track in the spectrometer (as- 
sociated to the proton) and two clusters in the electromagnetic calorimeter (associated to the 
photons from 7r° — > 77 decay) within a time window of 2 ns. The longitudinal position of the 
S + decay vertex was determined using the 7r° mass constraint to calculate the distance of its 
decay point from the calorimeter: 

A,,. = (1) 

where E± and E2 are the measured energies of the two clusters and r\i is the distance between 
the two clusters in the transverse plane. Good candidates were kept if the reconstructed p7r° 
invariant mass was within 6 MeV/c 2 of the nominal E + mass value. The mass interval was 
tightened from 8 MeV/ c 2 to 6 MeV/ c 2 with respect to the normalization channel (see below) to 
reduce the higher background contamination in the muon channel. 

Muon identification was achieved by requiring the presence of in-time signals from the first 
two planes of the muon detector (±2 ns with respect to the time measured in the charged ho- 
doscope). In addition, to reject pions and electrons, the energy deposited in the electromagnetic 
calorimeter in association to the muon track was required to be less than 2.5 GeV. 

The lower momentum threshold for the muon track was set to 7 GeV/c (it was 4 GeV/c 
for the electron channel) to reduce the background contamination and to increase the efficiency 
for muon reconstruction (see section [7|) . 

The muon momentum calculated in the S + rest frame was required to be less than 
0.125 GeV/c, exploiting the fact that no contribution is expected from the signal sample above 
this limit. This cut was not applied in the normalization channel. Similarly, since the proton 
momentum in the signal sample is mostly above 54 GeV/c, this criterion was used to enhance 
the probability that sufficient energy is deposited in the electromagnetic and hadron calorimeters 
to satisfy the trigger condition -E/HAC+LKr > 30 GeV. In the normalization channel the lower cut 
on the proton momentum was set at 40 GeV/c. 

The H° decay vertex position was obtained by computing the closest distance of approach 
between the extrapolated S + line-of-flight and the muon track. This distance was required to be 
less than 4 cm. Furthermore, the deviation of the transverse H° vertex position from the nominal 
line-of-flight defined by a straight line going from the center of the K$ target to the center of 
the liquid krypton calorimeter was required to be less than 3 cm. 

The longitudinal position of the H° vertex was required to be at least 6.5 m downstream 
of the Kg target, i.e. 0.5 m after the end of the final collimator and at most 40 m from the 
target. Similarly, the S + vertex position was required to be at least 6.5 m downstream of the 
target but at most 50 m from the target. The latter value was chosen larger than the upper 
limit for the S° vertex position to account for the lifetime of the S + particle. The longitudinal 
separation between the H° and E + decay vertices was required to be between —8 m and 40 m. 
The negative lower limit, tuned with Monte Carlo events, was chosen such as to take properly 
into account resolution effects. 

The quantity r*coG was defined as r*coG = Yli^i-^i/Yli^i where Ei is the energy of 
the detected particle and r*j the corresponding transverse position vector at the liquid krypton 
calorimeter position z^xr- For a charged particle, the quantity r*j was obtained from the extrap- 
olation to zlkt of the upstream segment of the associated track. For kinematical reasons, the 
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missing transverse momentum (pt) is smaller in the muon case with respect to the electron case. 
Therefore fcoG was required to be less than 8 cm instead of 15 cm as for the electron channel. 

By requiring the invariant mass ir + ir fi~ to be less than 0.490 GeV/c 2 , the contamination 
from Kl — > 7r + 7r~7r°, when the tt~ is misidentified as a muon, was reduced to a negligible level. 
This cut was not applied in the normalization channel. 

Cuts were also applied on the positions of the hit points of the tracks in the chambers 
and on the cluster positions in the electromagnetic calorimeter to improve the trigger and the 
reconstruction efficiencies. Furthermore the energies of the photons coming from the 7r° decay 
were requested to be between 3 and 100 GeV to ensure linearity on the LKr measurement. 

With the above selection criteria, 97 H° — > S + // - I7„ candidates were observed in the signal 
region. The distribution of events in the pir° invariant mass variable is shown in Figure Q] after 
all selection cuts were applied. Signal events peaking around the S + mass are clearly visible 
above the background. 

A contribution to the background (about 20% of the total) comes from overlapping events 
in the detector (accidentals). This contribution was estimated directly from data samples, looking 
to the activity in the detectors not in time with the main event time. There is a small contribution 
to the background from the decay H° — > Att° (populating the left side of the p — 7r° invariant 
mass distribution) with A — > pir~ and with ir~ either misidentified as a muon or decaying into 
7r~ —7- pTVp,. This contribution was estimated by the Monte Carlo simulation. However the main 
contribution to the background is due to scattered events in the final collimator of the neutral 
beam, in analogy to what was seen in the normalization channel. 




Invariant p%° mass (GeV/c 2 ) 

Figure 1: Reconstructed p7r° invariant mass distribution for S° — > £ + [x~V^ candidates after 
all selection criteria were applied. Points with error bars are data. The peak at the S + mass 
value shows clear evidence for the signal. The vertical dashed lines delimit the signal region. The 
background was evaluated with a likelihood fit of the data performed in two intervals, between 
1.164 and 1.180 GeV/c 2 and between 1.198 and 1.240 GeV/c 2 . The solid histo gram shows the 
sum of the background contribution (evaluated from the fit) and the Monte Carlo sample of the 
signal normalized to the events found in the data after background subtraction. 

Due to the difficulty to simulate this contribution and due to the low statistics in the 
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control samples coming from data, the background distribution in the pir invariant mass was 
fitted with an exponential in the intervals 1.164 - 1.180 GeV/c 2 and 1.198 - 1.240 GeV/c 2 . The 
fit was giving an estimate of (30.8 ± 3.8 s t a t ± 1.9 syst ) background events when extrapolated into 
the signal region. The systematic uncertainty on the background level was estimated by varying 
the fitting function and the fit region (also including the signal region, with the signal fitted 
with a Gaussian distribution). 

The data sample for the normalization channel H° — > T l + e~u e consists of 6316 events with 
a background of (3.4 ± 0.7)%. 

A detailed description of the reconstruction and selection for the normalization channel 
is reported in [7]. For that decay, since the electron is completely absorbed in the LKr, the 
corresponding track was identified by requiring a ratio between the energy deposit in the LKr 
and the momentum measured by the spectrometer (E/p) greater than 0.85 and lower than 
1.15. The other differences in the selection criteria of the signal and normalization channels are 
described above. 

6 Acceptance 

The acceptance for both signal and normalization decay channels was computed using a 
detailed Monte Carlo program based on GEANT3 [9l[T0]. Particle interactions in the detector 
material as well as the response functions of the different detector elements were taken into 
account in the simulation. A detailed description of the generator of the electron channel can 
be found in [7j. The generator for the muon channel was modified to include the contribution 
from pseudo-scalar currents [11] . parameterized with the form factor g% which, under Partially 
Conserved Axial Current (PCAC) hypothesis, can be extracted at q = from the Goldberg- 
Treiman relation [X2^ [T3] : 

S8(0)//i(0) = 2(M EO /M K -) 2 g l (0)/f l (0). (2) 

Since the g$ term is multiplied by mi e pton/^s ; its contribution is non-negligible for the muon 
case [llj . Using the available experimental results [7] 114} [T5] for the electron channel, the best 
estimates for the remaining non-vanishing form factors are: 

f2(q 2 = 0)/fi(q 2 = 0) = 2.0 ±1.3 

9i(q 2 = 0)/h(q 2 = 0) = 1.21 ± 0.05. (3) 

The central values were plugged into the Monte Carlo generator and the corresponding errors 
were used to evaluate the systematic error related to the acceptance calculation. Radiative 
corrections were not included in the generator of the muon channel. This leads to a systematic 
uncertainty of 1%, estimated using the Monte Carlo simulation for the electron channel with 
the electron mass substituted by the muon one. The acceptance for the signal H° — > S + 
was calculated to be (3.17 ± 0.01)%, while the acceptance for the normalization H° — > Y, + e~V e 
was (2.49 ± 0.01)%. Both quoted uncertainties originate from the statistics of the Monte Carlo 
samples. 

7 H° — > Y< + [i branching ratio 

The H° — > Y, + n^V^ branching ratio was obtained from the background-subtracted num- 
bers of selected events for signal and normalization, the corresponding acceptance values, the 
normalization branching ratio [7j and the efficiency on muon identification. These quantities are 
summarized in Tableland yield: 

BR(~° -> = (2-17 ± 0.32 stat ± 0.17 syst ) x 10~ 6 , (4) 

where the statistical uncertainty originates from the event statistics and the systematic one is 
the sum in quadrature of the various contributions presented in Table [2j 
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Table 1: Parameters used for the BR(H° — > 
normalization channel are taken from reference W 



measurement. The numbers used for the 









Event statistics 
Background 
Acceptance 
Muon inefficiency 
Branching ratio 


97 

(30.8 ± 4.2) events 
(3.17 ±0.01)% 
(1.5 ±0.5)% 


6316 
(3.4 ± 0.7)% 
(2.49 ±0.01)% 

(2.51 ± 0.03 stat ± 0.09 syst ) x 10" 4 



The largest contribution to the total systematic uncertainty comes from the background 
subtraction, described above. 

From the systematic uncertainty related to the measurement of the branching ratio of the 
normalization channel, the trigger efficiency contribution was eliminated, since it is common to 
both channels. A further systematic of 3% was added to take into account the dependence of 
the trigger efficiency on the lepton momentum. 

The sensitivity of the branching ratio measurement to the form factors was studied by vary- 
ing <7i(0)//i(0) and /2(0)//i(0) within the limits provided by their uncertainties and doubling or 
neglecting the 53(0) value. The muon momentum distribution from Monte Carlo simulation was 




Muon momentum 



(GeV/c) 5 



Figure 2: Efficiency for muon identification as a function of muon momentum (right), the dashed 
line shows the lower threshold at 7 GeV/ c applied to muon momentum. 

divided by the distribution of muon efficiency as a function of muon momentum as measured 
from — > [i^Vp (Kfj,2) decays, obtained from a large sample collected in 2003 [16] (see Figure 
[2]). A consistent result was obtained by considering K° — > ir^fi^u decays collected in 2002 but 
with much lower statistics. The resulting correction of (+1.5 ± 0.5)% was applied in the BR 
calculation. 



8 Conclusion 

Using data collected in 2002 with the NA48 detector at CERN, we obtain clear evidence 
of the decay H° — > 



Ua, with a precision on the branching ratio being significantly better 
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Table 2: Sources of systematic uncertainties. 



bource 


Uncertainty 


Background 


±0.4% 


AT 1 " J. " 

JNormahzation 


±3.0% 


L2 trigger efficiency 


±o.U/o 


Form factors 


-1-1 


Radiative corrections 


±1.0% 


Muon reconstruction efficiency 


±0.5% 


Total 


±7.9% 



than the existing published value: 

BR(~° -> E+zx-P^) = (2.17 ± 0.32 stat ± 0.17 syst ) x 10~ 6 . (5) 

This result is in good agreement with the branching ratio measured by the NA48 / 1 collaboration 
for the electron channel, once the theoretical ratio of the corresponding decay amplitudes is taken 
into account [3]. 
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